Submitted to ApJ: 9 November 1999; Revised: 24 January 2000 



Stellar Dynamics and the implications on the merger evolution in NGC 6240 

M. Tecza, R. Genzel, L. J. Tacconi, S. Anders, L. E. Tacconi-Garman and N. Thatte 
M ax- Planck- Institut fur extraterrestrische Physik (MPE) 
D-85748 Garching, Germany 

ABSTRACT 

We report near-infrared integral field spectroscopy of the luminous merging galaxy 
NGC 6240. Stellar velocities show that the two K-band peaks separated by 1'.'6 are the 
central parts of inclined, rotating disk galaxies with equal mass bulges. The dynamical 
masses of the nuclei are much larger than the stellar mass derived from the K-band 
light, implying that the progenitor galaxies were galaxies with massive bulges. The 
K-band light is dominated by red supergiants formed in the two nuclei in starbursts, 
triggered « 2 x 10 years ago, possibly by the most recent perigalactic approach. Strong 
feedback effects of a superwind and supernovae are responsible for a short duration 
burst (~ 5 x 10 6 years) which is already decaying. The two galaxies form a prograde- 
retrograde rotating system and from the stellar velocity field it seems that one of the 
two interacting galaxies is subject to a prograde encounter. Between the stellar nuclei is 
a prominent peak of molecular gas (H2, CO). The stellar velocity dispersion peaks there 
indicating that the gas has formed a local, self-gravitating concentration decoupled from 
the stellar gravitational potential. NGC 6240 has previously been reported to fit the 
paradigm of an elliptical galaxy formed through the merger of two galaxies. This was 
based on the near-infrared light distribution which follows a r 1//4 -law. Our data cast 
strong doubt on this conclusion: the system is by far not relaxed, rotation plays an 
important role, as does self-gravitating gas, and the near-infrared light is dominated by 
young stars. 

Subject headings: galaxies: individual (NGC 6240) — galaxies: interactions — galaxies: 
kinematics and dynamics — galaxies: starburst — infrared: galaxies 

1. Introduction 

The luminous infrared galaxy (LIRG) NGC 6240 (Wright, Joseph & Meikle 1984; Thronson 
et al. 1990) has a remarkable disturbed morphology which is characterised in the visible by loops, 
branches and arms extending out to 50kpc (Fosbury & Wall 1979). From this large-scale mor- 
phology, and the discovery of two nuclei in the central region Fried & Schulz (1983) concluded 
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that NGC 6240 is an interacting and merging system of two galaxies. This conclusion is sup- 
ported by spectroscopic observations in the visible (Fosbury & Wall 1979; Fried & Ulrich 1985) and 
near-infrared wavelength range (Herbst et al. 1990; Van der Werf et al. 1993; Sugai et al. 1997). 
NGC 6240 is the most luminous source of near-infrared line emission (L(H2)tot ~ 1O 9 L , Wright 
et al. (1984)). The H 2 emission is most likely excited in shocks triggered by the collision of the two 
galaxies. 

NGC 6240 has a bolometric luminosity L ho \ « L m = 6 x 10 11 L (for H = 75kms~ 1 Mpc -1 , 
D = 97 Mpc). This bolometric luminosity classifies NGC 6240 almost as an ultra luminous infrared 
galaxy (ULIRG; Lir > 1O 12 L ). LIRGs and ULIRGs emit more than 90% of their bolometric 
luminosity in the infrared wavelength range. There is consensus that the infrared emission comes 
from warm dust but what is the primary energy source responsible for heating the dust is still a 
matter of debate. The two possibilities are a deeply dust embedded active galactic nucleus (AGN) 
or a super starburst. Observations with the Infrared Space Observatory (ISO) (Genzel et al. 1998) 
show that the majority of ULIRGs, including NGC 6240, are starburst dominated. On the other 
hand, detection of a highly absorbed hard X-ray source indicates that NGC 6240 also contains a 
powerful AGN (Vignati et al. 1999). Most ULIRGs are interacting or merging systems with the 
star formation and/or AGN activity triggered by gas compression and inflows towards the center 
(Sanders et al. 1988). As such NGC 6240 is a nearby test case of a luminous merger. It can 
be used to test predictions of interacting and merging galaxies (Barnes & Hernquist 1996; Mihos 
& Hernquist 1996). The effect of the interaction on the star formation can be compared with 
model results to assess the true nature of ULIRGs. To investigate the stellar content and the 
properties of the starburst of NGC 6240, to determine the stellar dynamics and to test predictions 
of interacting galaxy models we have carried out high resolution near-infrared spectroscopy with 
the Max-Planck-Institut fur extraterrestrische Physik (MPE) integral field spectrometer 3D. 



2. Observations and Data Reduction 

NGC 6240 was observed with the MPE near-infrared imaging spectrograph 3D (Weitzel et 
al. 1996) in conjunction with the tip-tilt correction adaptive optics system ROGUE (Thatte et 
al. 1995) in two observing runs. 3D is an integral field spectrograph that simultaneously obtains 
spectra for each of 256 spatial pixels covering a square field of view with over 95% fill factor. In 
both observing runs the spectral resolving power (R = A/ A A) was 2000 and Nyquist sampled using 
two settings of a piezo-driven flat mirror. 

The first observing run took place in April 1996 at the ESO 2.2 m telescope on La Silla, Chile. 
The pixel scale was tf!3 per pixel, and a wavelength range from 2.18/xm to 2.45^m was covered. 
The total on-source integration time was 4600 s with individual frame integration times of 300 s or 
400 s. The same amount of time was spent off-source 1' E and W of the nuclear region of NGC 6240 
for sky background subtraction. The seeing during the observations was better than 0'/8. 
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The second observing run took place in March 1998 at the Anglo Australian Telescope in 
Coonabarabaran, Australia. Here the pixel scale was (X'4 per pixel and a wavelength range from 
2.155/um to 2.42^im was covered. This wavelength setting includes the ~H.2V = 1 — S(l) emission 
line at a redshifted wavelength of 2.173/im. The total on-source integration time was 3100 s with 
individual integrations of 100 s each. Again, the same amount of time was spent off-source 1' E 
of the nuclear region of NGC 6240 to subtract the sky background. The seeing throughout the 
observations was better than 1". 

To correct for the atmospheric transmission a reference star was observed before and after the 
science integrations. The data were reduced using the set of 3D data analysis routines written 
for the GIPSY (Van der Hulst et al. 1992) data reduction package. This included wavelength 
calibration, spectral and spatial flat fielding, dead and hot pixel correction, and division by the 
reference stellar spectrum. Data cubes from the individual exposures were recentered and added 
using the centroid of the broadband emission from the southern nucleus. Absolute flux calibration 
was established by comparison of the broadband emission with K-band photometry by Thronson 
et al. (1990). Emission line and absorption line maps as well as the continuum map were extracted 
by performing a linear fit through line-free regions in the vicinity of the lines in the spectrum of 
each spatial pixel. 



3. Results 

3.1. Line Maps and Spectra 

The spectrum of NGC 6240 in the range from 2.15/xm to 2.45/um contains a number of molecular 
and atomic emission lines, as well as stellar absorption features. 

Figure 1 shows the line-free continuum map in the wavelength range from 2.2//m to 2.45/um, 
a Brackett 7 map and a map of the CO 2 — absorption bandhead. All three maps are dominated 
by two compact sources (hereafter referred to as nuclei) separated by l'/6. The double nuclear 
structure has been observed previously in the visible (B, V, R) (Keel 1990) and NIR (J, H, K) 
bands (Thronson et al. 1990) but the maps are also similar to what is observed at other wavelengths 
(Condon et al. 1982; Colbert, Wilson & Bland-Hawthorn 1994; Keto et al. 1997). In addition to the 
double nucleus the Brackett 7 map shows emission extended out to l'/5 northwest of the northern 
nucleus in the direction of the third radio emission peak N3 observed by Colbert et al. (1994). 

Also shown in Figure 1 is the integrated H2t> = l — S(l) emission. In contrast to the K- 
band continuum the H2 v = 1 — 0S(1) emission shows a single emission peak between the nuclei with 
filaments extending from it. The emission peak is at a separation of C//35 from the southern nucleus 
towards the northern nucleus and is extended in the northern direction. This H2 morphology has 
been previously observed with Fabry- Perot imaging techniques by Herbst et al. (1990), Van der 
Werf et al. (1993) and Sugai et al. (1997). The H2 morphology is very similar to the CO J = 2 — 1 
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morphology, indicating that also the near-infrared H2 emission lines trace the bulk of the molecular 
material in NGC 6240. 

The line maps in Figure 1 show three major regions of interest. These are the two nuclei, and 
the region of the H2 v = 1 — S(l) emission peak between them. Figure 2 shows the spectra of these 
regions, each within an 1" diameter circular aperture. Even though the peak of the H2 emission 
lies between the nuclei, H2 dominates the line emission in all three regions. In fact, NGC 6240 
has the highest ratio of H2 t> = l — S(l) luminosity to bolometric luminosity of any galaxy which 
has thus far been observed (Wright et al. 1984). The H2 lines are very broad and have a full 
width half maximum (FWHM) of 550 km s -1 with line wings extending over 1600kms _1 full width 
zero power (FWZP). The emission line profiles are asymmetric with red and/or blue wings whose 
strength varies over the emission region. Figure 3 shows a comparison of the line-of-sight velocity 
profiles of the H2 v = l — S(l) and CO J = 2 — 1 lines for the nuclear region of NGC 6240. 

A characteristic of NGC 6240 is its low Brackett 7 equivalent width (e.g. Lester, Harvey & 
Carr 1988; Van der Werf et al. 1993). In earlier work (see references in Van der Werf et al. 1993) 
the line was not or only marginally detected. In the improved sensitivity 3D data the line is now 
clearly detected. The Brackett 7 equivalent width of the northern and southern nucleus is 6.5±1.6A 
and 2.9 ± 0.4A, respectively. A line map and spectra in different positions can be extracted from 
the data cube. The FWHM of the Brackett 7 line is 716 km s -1 and 689 km s -1 for the northern 
and southern nucleus, respectively, measured over a 1" diameter circular aperture. The northern 
nucleus is redshifted with respect to the southern nucleus by 137 ± 49kms _1 , a value similar to 
the velocity difference of 147kms _1 measured by Fried & Ulrich (1985) with optical emission and 
absorption lines. 

The 3D spectra also show that NGC 6240 has deep CO absorption bandheads with a large 
velocity dispersion. From similar data Lester & Gaffney (1994) and Doyon et al. (1994) deduce a 
velocity dispersion of 350 km s _1 and 359 km s _1 , respectively, for the region of NGC 6240 containing 
both nuclei. Using 3D, we are able to resolve spatially the velocity dispersions of each nucleus. 
A detailed analysis of this is presented in §5. The equivalent width of the CO 2 — bandhead is 
determined using the wavelength intervals given in Origlia, Moorwood & Oliva (1993) and corrected 
for the velocity broadening using the formula in Oliva et al. (1995). The velocity dispersions arc 
the ones derived in §5 and the resulting equivalent widths are 13 ± 3A and 15 ± 1 A for the northern 
and southern nucleus, respectively. 



3.2. Extinction 

From the spectra in Figure 2 it is apparent that the spectrum of the H2 emission peak has 
a shallower, i.e. redder, slope than the two nuclei. Assuming an intrinsically constant continuum 
slope over the entire nuclear region of NGC 6240, this reddening can either be due to extinction 
or a nonstellar contribution to the spectrum at long wavelengths. The spectral energy distribution 
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(SED) of NGC 6240 (Draine & Woods 1990) shows no nonstellar continuum below 5/im. Above 
5/iin thermal emission from dust with a temperature of ~ 200 K is dominant. Different spectral 
slopes over the nuclear region of NGC 6240 thus are almost certainly due to extinction variations. 

To derive the extinction values the spectral slope of all pixels in the data cube are compared 
with the spectral slope of a K4.5 supergiant, which was observed with 3D during the La Silla 1996 
observing run. A direct comparison of stellar spectra with the spectrum of NGC 6240 yields a K4.5 
supergiant as the best fitting template (see §4.1). With the interstellar extinction law from Draine 
(1989) 

^ = 0-351 • 

and an uniform foreground screen model (UFS) we derive the extinction map shown in Figure 4. 
The extinction towards the southern and northern nucleus is Ay = 5.8 and Ay = 1.6, respectively. 
The peak extinction is Ay = 7.2 ± 0.7. For a mixed model, where absorbing dust and stellar 
emission are completely mixed, the peak value is Ay = 18.4^'g. The morphology of the extinction 
map with a single peak between the continuum nuclei indicates a dust concentration between the 
nuclei and is coincident with the H2^ = l — 0S(1) and CO J = 2 — 1 emission peaks. 



4. The Infrared Nuclei 

4.1. The Nature of the K-band Continuum 

To determine the stellar type dominating the near-infrared light in starburst galaxies the 
CO 2 — absorption bandhead is an often used spectral absorption feature. For late- type giants 
and supergiants this absorption feature is rather deep, and even with data of medium spectral 
resolution the equivalent width can be used to determine the stellar type. There is an ambiguity in 
the CO 2 — equivalent width between KI and Mill stars, however. From the CO 2 — absorption 
bandhead equivalent width the spectrum of NGC 6240 can be explained either by a population 
of red giants or red supergiants. However, because of the high spectral resolution of our data, 
we can resolve the discrepancy by a direct comparison of the spectrum of NGC 6240 with both 
stellar types. Figure 5 compares the spectrum of the CO absorption bandheads of the southern 
nucleus of NGC 6240 with four stellar spectra of giants and supergiants also observed with 3D at 
the same spectral resolution (Schreiber 1999). For the comparison the spectrum of NGC 6240 was 
shifted into the restframe of zero redshift, and the stellar spectra were broadened with a Gaussian 
velocity profile of 600 km s -1 FWHM. Clearly the K and M giant spectra do not fit the deep 
absorption features and the continuum slope. On the other hand the absorption bandheads of a 
M3+ supergiant are too deep to represent the NGC 6240 spectrum. The best matching stellar 
spectrum is that of a K4.5 supergiant. However, its CO absorption bandheads are still somewhat 
too shallow and we conclude that the near-infrared light is due to late K (K5) or early M (M0/M1) 
supergiants. Our conclusions are consistent with those of Sugai et al. (1997) whose data have lower 
spectral resolution but cover a wider spectral range. The fact that the K-band light of the nuclei 
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is dominated by red supergiants indicates a starburst as the source of the K-band luminosity. The 
presence of supergiants implys that the starburst was triggered quite recently. 



4.2. Starburst Simulations 

To further constrain the age and other characteristics of the starburst we used the program 
STARS (Schreiber 1999) to simulate several starburst scenarios. In Figure 6 we show model pre- 
dictions for the variations of Brackett 7 equivalent width, the K-band luminosity to stellar mass 
ratio (Lk/M*) and the starburst bolometric luminosity to K-band luminosity ratio (L£ o1 /Lk) for 
four burst durations: 1, 5 and 20 million years as well as continuous star formation. We adopted 
a Salpeter (1955) initial mass function (IMF) for masses between 100 M® and 1M®. The calcu- 
lated luminosity output of the starburst is insensitive to the shape of the IMF at the low mass end 
(< 1 A4®). The low mass stars do contribute to the total mass of the stars formed in the starburst, 
however. To take this into account, we integrated a Miller-Scalo IMF (Miller & Scalo 1979) from 
IA4® down to 0.08 M®, the lower mass limit for hydrogen burning. Because stars with masses 
> 25 M.® do not evolve into red supergiants, the K-band luminosity does not depend on the choice 
of the upper mass cutoff, and due to their small number, their contribution to the total mass is 
negligible. 



4.3. Starburst Age and Duration 

The CO 2 — absorption bandhead equivalent width (or better the presence of late type su- 
pergiants) allows us to determine the age of the starbursts in NGC 6240. The late K or early M 
supergiants that dominate the K-band light of the nuclei have masses from 10 to 20 A4© and a 
typical age of 15 to 25 million years. Since red supergiants dominate the K-band luminosity only 
during this period, the starburst must be of similar age. This age is shown in Figure 6 as a vertical 
hatched bar. 

The Brackett 7 emission line equivalent width allows us to constrain the duration of the 
starburst activity. The agreement of the Brackett 7 and K-band continuum morphologies indicates 
that both the line and continuum emission originates in the same star forming regions. The H2 
emission can be explained by excitation in slow C-shocks (Van der Werf et al. 1993; Sugai et al. 
1997; Egami 1998). Since the morphology of the H2 emission differs from that of the Brackett 7 
emission, shocks most likely do not contribute to the Brackett 7 emission. In a starburst, only 
stars with masses > 20 — 30 M.® ionise their surrounding medium to produce Brackett 7 emission, 
but due to their high mass their lifetime is very short (< 10 million years). The low Brackett 7 
equivalent width in NGC 6240 indicates a low number of hot, young stars, either because the 
starburst is aging or because stars with masses > 20 M® were never formed. We prefer the former 
explanation, because of the signatures for 50 to 100 M.® stars in nearby starburst templates (see 
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Thornley et al. 1999). 

The middle panel of Figure 6 shows that the starburst duration is < 5 million years. The short 
starburst duration can be explained by strong negative feedback effects from the starburst itself. 
The onset of vigorous star formation produces young, massive stars and, thus, supernovae in a very 
short time (< 5 million years). The winds of these stars deplete the molecular gas and the star 
formation subsides. The existence of a superwind in NGC 6240 was shown by Heckman, Armus &; 
Miley (1990) from Ha line mapping and spectroscopy. 

The duration of the starbursts is much smaller than the age of the starburst but comparable 
to the dynamical time scale of ~ 7 million years of the two rotating nuclei (see §5.1). The time 
scale for the interaction and merging of the two galaxies is several hundred million years and hence 
much larger than the starburst duration and starburst age. Such short but violent star formation 
events are predicted by models of interacting and merging galaxies, where they are triggered by 
the close encounters of the interacting partners (Mihos & Hernquist 1996). 

4.4. Stellar Mass in the Starburst 

The light-to-mass ratio (Lk/M*) of the starburst can be calculated from simulations of the 
K-band luminosity and the total mass of stars formed in the starburst. The middle panel of Figure 
6 shows Lk/M* as a function of starburst age. It peaks at the time when the red supergiants 
dominate the K-band luminosity. For the starburst age in NGC 6240 the simulation yields 1 < 
Lk/M* < 3. From the dereddened K-band luminosities in Table 1 within 1" a stellar mass 
M* = 0.4 — 1.2 x 10 8 A4© and A4* = 0.8 — 2.3 x 10 8 Mq is derived for the northern and southern 
nucleus, respectively. Here the lower mass cutoff was assumed to be 0.08 Mq. If the lower mass 
cutoff was 0.25 Mq, 0.5 Mq or IMq the estimated stellar mass would be smaller by 12%, 23% 
or 39%, respectively. This is the mass of young stars formed in the starburst contained within the 
central 500 pc of the two nuclei. It shows that the two nuclei are massive objects and possibly could 
be the nuclei of two galaxies involved in a collision, based on the mass of young stars alone. 

4.5. The Luminosity of the Nuclei 

The bottom panel in Figure 6 shows the starburst bolometric luminosity to K-band luminosity 
ratio (L£ o1 /Lk) as a function of starburst age. For a starburst age of 15 to 25 million years this ratio 
is ~ 100. From the dereddened K-band luminosity within 2" diameter of each nucleus (see Table 
1) this ratio yields L£ ol s=s 0.7 x 10 11 L Q and L£ ol « 1.7 x 10 11 L Q for the northern and southern 
nucleus, respectively, totaling to L£ ol « 2.4 x 1O 11 L . This is | to \ of the entire bolometric 
luminosity L^ Q \ «6x 10 11 L Q of NGC 6240 in the IRAS beam. Within a 5" diameter aperture 
the total K-band luminosity Lk = 6.1 x 10 9 L Q is ~jgo °^ tne total bolometric luminosity. In the 
bottom panel of Figure 6 this ratio is shown as a horizontal, hatched bar. Within the errors of the 
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simulations the entire bolometric luminosity of NGC 6240 can be explained by the starburst. 

The claim that the starbursts in the nuclei contribute significantly to the bolometric luminosity 
of NGC 6240 is further strengthened by radio observations. The radio flux and the far infrared 
(FIR) luminosity of NGC 6240 follow the radio-to-FIR correlation for starburst galaxies. Lisenfeld, 
Voelk k, Xu (1996) determined empirically the parameter 



relating the radio-power to the FIR-luminosity. For starburst galaxies q2.4GHz = 2.40±0.22. Taking 
the radio-power S2.4 GHz = 1-7 x 10 23 WHz -1 at v = 2.4 GHz for the nuclear region including both 
nuclei, and comparing it with the FIR-luminosity S Fm = 1.9 - 2.3 x 10 38 W yields for NGC 6240 
Q2.4 GHz = 2.47 — 2.55, within the uncertainties of the radio-to-FIR correlation. If we include the 
radio-emission from the western region (region W in Colbert et al. 1994), which has no counterpart 
in the visible or infrared, we get q2.4GHz = 2.31 — 2.39, also within the uncertainties of the q- 
parameter for starburst galaxies. 

X-ray observations and the detection of 25.9/mi [O IV] emission indicate the presence of a 
powerful AGN in NGC 6240 as well (Schulz et al. 1998; Komossa, Schulz & Greiner 1998; Iwasawa 
& Comastri 1998; Genzel et al. 1998; Vignati et al. 1999). In fact, Vignati et al. (1999) conclude from 
their BeppoSAX data that the highly absorbed (Nh ~ 2 x 10 24 cm~ 2 ) AGN can account for most 
of the bolometric luminosity of NGC 6240. These two conclusions seem to contradict each other 
. However, extrapolating from the near /mid-infrared spectroscopy and hard X-ray photometry to 
bolometric luminosity and the possibility of significant correction of the observed luminosity for 
anisotropy introduce uncertainties of factors > 2. We conclude that NGC 6240 is a composite 
object where both star formation and AGN play a role. 

At the position of the K-band nuclei also Brackett 7 emission is observed, however, very 
little continuum emission is detected at the position of the north-west extension of the Brackett 7 
emission. This could be an indication that the Brackett 7 emission is arising from a very young 
starburst that has not yet formed enough stars to contribute to the continuum emission, but whose 
hot, young stars have already ionised the interstellar medium (ISM). If true, we can estimate the 
bolometric luminosity of this emission region using the radio-to-FIR correlation and get S'fir = 
0.4 x 10 11 L Q . This is about half the bolometric luminosity of the northern nucleus and « 15% of 
the total bolometric luminosity of NGC 6240, a non-negligible contribution to the total bolometric 
luminosity of NGC 6240. Colbert et al. (1994) have interpreted their third radio peak N3, based on 
its steeper radio spectrum, as a clump of electrons, driven away from the nucleus by a superwind. 
From the Brackett 7 emission we think this explanation is less likely and favour instead that N3 is 
a young starburst region. 
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5. Stellar Kinematics 

In the near-infrared wavelength range from 2 — 2.5/xm stellar kinematics can be determined 
from absorption features of Ca, Na and CO. By far the strongest absorption feature is the CO 2 — 
absorption bandhead at 2.29/mi. It has a very sharp blue edge which is very sensitive to stellar 
motions. Another advantage is that the extinction in the K-band is only 10% of the visual extinc- 
tion. The CO 2 — absorption bandhead is, therefore, very well suited for the determination of the 
stellar kinematics in the gas- and dust-rich infrared galaxies. 

We have derived the velocity dispersions and radial velocities of the nuclei of NGC 6240 by 
using a modified Fourier correlation quotient (FCQ) method (Bender 1990; Anders 1999) and a 
direct fitting method on the CO 2 — absorption bandhead data. In the FCQ method the cross- 
correlation of the template spectrum with the galaxy spectrum and the auto-correlation of the 
template spectrum are computed in the Fourier domain. Then the quotient of the cross-correlation 
and the auto-correlation is calculated. To suppress high-frequency noise a Wiener filter is applied 
to the correlation-quotient. Fourier transformation of the correlation quotient finally yields the 
line-of-sight velocity distribution in the galaxy spectrum. In the direct fitting method a stellar 
template spectrum was convolved with a Gaussian broadening function. The broadening function 
was parameterised by a radial velocity and a velocity dispersion. With a least-x 2 fitting technique 
the best fitting set of parameters was computed. To derive the kinematic parameters with the FCQ 
method the deconvolved velocity profile is fit with a Gaussian profile. For the template spectra used 
in both methods we chose a K4.5 supergiant spectrum from the 3D stellar library (Schreiber 1999). 
While the FCQ method is rather insensitive to template mismatch, the results from the direct 
fitting method depend very much on the choice of the stellar template spectrum. We find that the 
best choice for the type of the stellar template is a late K or an early M supergiant (see also §4.1). 
The errors of both the FCQ and the direct fitting method were determined by modeling of galaxy 
spectra from stellar spectra with different noise levels and trying to recover the input spectrum. 
With the direct fitting method a x 2 e( r anar y s i s was use< ^ as an independent way to determine the 
fit error. 



5.1. Stellar Mean Velocity field: Rotation and Dynamical Mass 

Because of the sharp edged profile of the CO absorption bandhead we were able to determine 
the stellar velocity field for the entire central region of NGC 6240 with the direct fitting method. 
The resulting velocity field is shown in contour form in Figure 7 superposed on the NICMOS ri- 
band HST archival image (PI: N. Scoville, proposal ID: 7219). Velocity gradients across both nuclei 
are clearly recognizable. The velocity error at the position of the southern nucleus is ±18kms _1 
and increases to ±50kms -1 at a distance of 076; for the northern nucleus the value is ±60kms _1 
and ±85kms _1 at a distance of 075. Along the edge of the velocity field shown in Figure 7 the 
uncertainty is ±150 km s -1 . 
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The southern nucleus shows a velocity gradient at a position angle of 34° west of the north- 
south axis with the south-east of the nucleus being redshifted. The velocity gradient of the northern 
nucleus has a position angle of 41° east of the north-south axis. The north-east of the northern 
nucleus is redshifted with respect to the center. The maximum relative velocity shift from the 
center along the velocity gradient within 1" is ±165 km s -1 and ±170 km s" 1 for the southern and 
northern nucleus, respectively. The agreement of the K-band morphology from NICMOS with the 
kinematical morphology of the stellar velocity field is remarkable. The two nuclei are elongated 
along the velocity gradient in the velocity field, indicating individual inclined rotating disks in the 
nuclei. The northern nucleus is redshifted with respect to the southern nucleus by ~ 50kms _1 . 
This value is smaller than the value measured in the Brackett 7 line and measured by Fried & 
Ulrich (1985) but consistent within the uncertainties of our measurement. 

The stellar velocity is surprisingly different compared to the velocity field of the molecular gas. 
The gas forms a rotating disk between the two nuclei (Tacconi et al. 1999) with a axis of rotation 
not aligned with either of the stellar disks in the nuclei. This implys a decoupling of the stellar 
motions from the gas motions, most likely caused by the tidal forces of the interaction. 

The dynamical mass of the NGC 6240 nuclei can be determined from the stellar velocity field. 
Under the assumption that the stars move on circular orbits around the center of each nucleus, the 
dynamical mass within the radius R for a inclination corrected rotation velocity v rot is given by 



The measured recessional velocities for the two nuclei are shown in Figure 8. We fitted a model 
rotation curve to the data with a solid body (v oc R) part for radii R < 100 pc and a flat part 
at a rotation velocity v Iot for larger radii. The model disk was inclined towards the line of sight 
by an angle of 45° and 60° for the northern and southern nucleus, respectively. The inclination 
angles were determined by fitting ellipses to the isophotes of the NICMOS H-band image. The 
model disk was convolved with a Gaussian profile with 0'/8 FWHM, corresponding to the average 
observing conditions. For the assumed inclinations the best fit of the rotation curve is achieved with 
v TOt = 270 ± 90 km s^ 1 for the southern nucleus and t> rot = 360 ±195 km s^ 1 for the northern nucleus. 
Even for an inclination of 90° the rotation velocities are 180 ± 90kms _1 and 200 ± 155 km s -1 , 
respectively. Assuming an extreme inclination angle, the dynamical mass within 500 pc of the 
northern nucleus is 2.4 x 10 9 M®, and 1.9 x 10 9 M® for the southern nucleus. Such large masses 
can only be explained if the observed K-band nuclei are indeed true nuclei of the two merging 
progenitor galaxies. Given the likely inclinations indicated by the NICMOS images the true masses 
are likely even factors > 2 higher. Table 2 lists the derived dynamical masses of both nuclei for an 
inclination of 45° and 60° for the northern and southern nucleus, respectively. 

We also list in Table 2 the ratio of dynamical mass to stellar mass content of the starburst 
derived from the light-to-mass ratio. The dynamical mass exceeds the visible mass by about an 
order of magnitude. The missing mass cannot be explained by a contribution of molecular gas at 
the position of the nuclei. The molecular gas is concentrated between the nuclei with a relatively 




(1) 
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small mass at the position of the nuclei (Tacconi et al. 1999). From the millimeter CO J = 2 — 1 
line flux within an 1" diameter aperture on the nuclei an upper limit for the cold molecular gas 
mass in the nuclei can be determined. For the northern and southern nucleus this upper limit is 0.2 
and 0.25, respectively, of the gas mass within the central 2" of the CO-disk. Taking 2 x 10 9 Mq for 
the mass of the gas concentration (Tacconi et al. 1999) the contribution from the cold molecular 
gas to the total mass of the nuclei is 4 — 5 x 10 8 M®, much smaller than the missing mass. The 
contribution from a dark matter halo to the bulge masses cannot be determined from our data. In 
normal spiral galaxies this contribution is negligible. The only remaining likely source of missing 
mass is in the old stellar population of the progenitor galaxies themselves. The number of faint 
stars could add up to a significant mass, yet their luminosity would be negligible compared to the 
red supergiants of the recent starburst. We conclude that the two K-band peaks are indeed the 
central bulges of the progenitor galaxies and that these bulges are fairly massive, in accordance 
with the prediction of Mihos & Hernquist (1996) and Barnes & Hernquist (1996). 



5.2. Stellar Velocity Dispersion: Mass Concentration between the Nuclei 

While the sharp edge of the CO 2 — bandhead allows a redshift determination even for data 
with a low signal-to-noise ratio, we cannot determine the velocity dispersion from such data. We 
have, therefore, completed the analysis of the spatial variation of the velocity dispersion with 
integrated spectra for several apertures centered in the nuclear region of NGC 6240. To determine 
the true velocity dispersion we first subtracted the stellar rotation component as derived above. 
Within the selected apertures we corrected the data cube for the stellar rotational velocity field 
by shifting the spectra of all spatial pixels to the same restframe. Figure 9 shows the distribution 
of the velocity dispersion over the nuclear region of NGC 6240. We show the calculated velocity 
dispersion for four apertures: two centered on the nuclei; one centered between the nuclei; and one 
south of the southern nucleus. The error bars along the declination axis denote the width of the 
apertures over which the spectrum was integrated. The aperture width in right ascension is 0.8". 
Figure 9 shows that the velocity dispersion peaks between the two nuclei at a = 276 ± 51kms _1 , 
while the velocity dispersion at the northern and southern nucleus is a = 174 ± 54kms _1 and 
a = 236 ± 24kms~ 1 , respectively. Because this result is derived from the velocity field corrected 
data cube, the broad velocity profile cannot be due to the overlapping of two velocity profiles from 
components separated by scales larger than the spatial resolution of our observations, that is 350 
to 400 pc. 

Replacing v TOt with the velocity dispersion a, equation (1) can also be used to estimate the 
mass M = 2.3 x 10 2 • R pc ■ cr^ m -i Mq within the radius R. This is a very simplistic form of 
the Jeans equation (Binney &: Tremaine 1987) where only the isotropic part is considered and 
numerical factors from the space distribution of the stars and the spatial dependency of rotation 
and velocity dispersion are neglected (thus resulting likely in an under-estimate of the dynamical 
mass). Following this argument the peak of the velocity dispersion indicates a mass concentration 
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between the two nuclei. We get masses for the northern and southern nucleus of 0.8 ±0.5 x 10 9 Mq 
and 1.8 ±0.4 x 10 9 Mq, respectively. Within the uncertainties this is consistent with the dynamical 
masses derived above. The mass between the nuclei is 2.1 ± 0.8 x 1$ M.q. The value of the 
masses determined in this way are somewhat questionable, since the Jeans equation applies to 
relaxed systems, a requirement that is not fulfilled in NGC 6240. At the peak of the stellar velocity 
dispersion little stellar continuum is detected. The low continuum flux between the nuclei cannot 
be explained by extinction which, although it peaks between the nuclei, is too low to explain a 
hidden concentration of stars. On the other hand, Tacconi et al. (1999) find a concentration of 
cold molecular gas of ~ 2 x 10 9 Mq centered within 1 kpc between the nuclei. It is thus very likely 
that the peak of the stellar velocity dispersion is caused by this molecular gas concentration. Our 
extinction and H2 v = 1 — S(l) maps indicate also that there is a corresponding peak of dust and 
hot H 2 . 



6. Molecular Gas 

6.1. Shock Excited H2 Emission 

From an analysis of the near-infrared H2 line ratios Van der Werf (1996) and Sugai et al. 
(1997) conclude that H2 is thermally excited with an excitation temperature of ~ 2000 K. The 
four H2 lines we observed also fit this result. With the Infrared Space Observatory (ISO) pure 
rotational H2 lines were detected for the first time in NGC 6240 (Lutz et al. 1996; Egami 1998; 
Rigopoulou et al. 1999). The analysis of their line ratios yields an excitation temperature < 400 K. 
Van der Werf (1996), Sugai et al. (1997) and Egami (1998) all find a slow continuous shock (C- 
shock) as the best fitting model for the H2 excitation. Fast shocks would lead to a discontinuous 
or jump shock (J-shock). Due to the discontinuity J-shocks ionise the interstellar medium and 
hence emission from ionised elements is expected. The low Brackett 7 to H2 v = 1 — S(l) line ratio 
therefore indicates strongly that the shock is non-dissociative. The agreement of the Brackett 7 
and continuum morphology suggests the starburst as the source of the Brackett 7 emission. 

In the C-shock models the shock has a velocity < 40kms~ 1 , yet, the observed H2 line widths 
are Av w 550 km s -1 . The large line width are most likely a superposition of several narrower 
H2 lines with different radial velocities along the line-of-sight. If we consider the distribution of 
molecular gas in form of dense clouds within a less dense intercloud medium we can explain the 
slow shock speed. The collision between molecular clouds occurs under a variety of angles and the 
shock speed therefore is only the projected radial velocity difference. But even for a direct collision, 
the shock speed in the denser cloud is smaller than in the lower density intercloud medium. An 
originally fast shock propagating in the intercloud medium will propagate with a slower shock speed 
when it enters a medium of higher density like a molecular cloud. It is these secondary slow shocks 
in the dense clouds that we are observing in the near-infrared ro-vibrational H2 lines. 
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6.2. Cold and Hot Molecular Gas 

The gas motions in the central region of NGC 6240 are highly turbulent. Tacconi et al. (1999) 
find CO J = 2 - 1 line width of up to 400 km s" 1 FWHM (and 1000 km s" 1 FWZP) and their channel 
maps show filaments extending from a central CO-disk out to 2kpc. 

A comparison of Figure 10 with the channel maps of the CO J = 2 — 1 line (Tacconi et al. 1999, 
Figure 2) indicates that the cold molecular gas (CO J = 2 — 1) and the hot, shock-heated H2 have 
a very similar overall morphology. Local differences between the CO J = 2 — 1 and H2 v = l — S(l) 
channel maps are apparent. The most obvious difference is in the central region between the nuclei. 
The hot gas is more extended than the cold gas and a disk structure as seen in the CO J = 2 — 1 
line is not visible. In addition, the H.2V = 1— 0S(1) line is 150 kms -1 to 250 kms -1 broader than 
the CO J = 2 — 1 line. Both effects can be explained with a different distribution of the cold and 
hot gas. Roughly half of the cold gas has already settled in the center between the nuclei, while 
the shock-heated interstellar medium of the two galaxies is still more extended. The H2 emission is 
excited on the surfaces of molecular clouds while the CO J = 2 — 1 emission originates more in the 
volume of the molecular cloud. Because we sample with the H2 emission a larger volume than with 
the CO J = 2 — 1 emission, we also sample a larger range in radial velocities and hence observe larger 
linewidths for the H2 lines. The larger line width of the H2 v = l — S(l) lines are not symmetric 
about the peak of the CO J = 2 — 1 lines but are blueshifted by 150 km s -1 . 

In the structure outside the nuclear region of NGC 6240 the H2 v = l — S ( 1 ) and C O J = 2 — 1 
emission is dominated by a filamentary structure. The filaments in the south-east and south-west 
are the most prominent and show similar morphology and kinematics in the CO J = 2—1 and 
H2 v = l — S(l) line. Also the weaker filaments north of the nuclei show this similarity. They are 
most likely gas flows towards the centers of the nuclei. 

A comparison of the velocity profiles of the CO J = 2 — 1 and H2 v = l — 0S(1) lines is shown 
in Figure 3. From the line-of-sight velocity profiles in Figure 3 several kinematical components 
of the H2 v = l — S(l) line are visible. They vary over the field, especially at the positions where 
the filaments extend from the H2 v = l — S(l) emission peak. North-east of the northern nucleus 
two components are recognizable with a velocity difference of s=s 250 km s" 1 . The filament in the 
south-west shows also two major components with a velocity difference of « 200 km s -1 . The 
H2V = 1 — 0S(1) line exhibits this multi-component character over the entire emission region of 
NGC 6240. 

The radial velocity difference between northern and southern nucleus is ps 150 km s^ 1 in the 
H2 v = l — S(l) line and « 100 kms -1 in the CO J = 2 — 1 line. This is similar to the radial velocity 
difference of the nuclei measured in the Brackett 7 line. Globally the molecular gas appears to 
follow the motion of the two nuclei. As with the cold gas, a rotation of the hot gas seems the 
best explanation for the velocity gradient of the H2 v = l — S(l) line. The velocity gradient peaks 
between the nuclei at the position of the CO J = 2 — 1 emission peak. At this position the velocity 
dispersion of the H2 v = l — S(l) emission also reaches its maximum of 240 km s -1 . 
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The central disk of cold molecular gas has a rotation axis which is not correlated with the 
rotation axis of either nucleus. Rather the sense of the rotation corresponds to the relative motions 
of the nuclei. The CO-disk seems to retain the memory of the orbital history of the interaction and 
is independent of the gravitational forces exerted by the nuclei. 



7. Interaction and Merging 

We have shown that NGC 6240 is a merging system of two galaxies with massive bulges. In 
the following sections results from our observations and published results are used to draw a more 
detailed picture of the interaction in NGC 6240. 



7.1. Prograde Encounter? 

The rotation velocities for both nuclei, derived in §5.1, are much larger than the relative 
velocity difference between the two nuclei. Assuming circular orbits for the nuclei the small radial 
velocity difference and projected distance between the nuclei yields a small projected orbital angular 
momentum of the nuclei. This means that the two galaxies either have lost a major fraction of their 
original orbital angular momentum in the interaction, or that the true orbital angular momentum 
is much larger than the projected value. A larger true orbital angular momentum means that the 
actual distance between the two nuclei is larger than the measured, projected separation. On the 
other hand, during the interaction of two galaxies their orbits can vary and the nuclei might rather 
be on radial than on circular orbits. In this case their orbital momentum could be still large. The 
exact trajectories of the galaxies depend on the initial conditions of the interaction. 

Because we can only observe projected angular momenta, the exact orientation of the disks 
with respect to the orbital plane cannot be determined. A schematic view of the collision geometry 
is shown in Figure 11. The projected spin angular momentum of the northern nucleus is almost 
parallel to the systems projected orbital angular momentum, while the projected spin angular 
momentum of the southern nucleus is nearly normal to it. If the projected angular momenta are 
the true angular momenta, one partner in the collision, now seen as the northern nucleus, was 
subject to a prograde encounter, while the other partner is inclined with respect to the orbital 
plane. For the motion of the two nuclei around each other we assume a circular orbit, whose 
position angle and inclination is equal to the position angle and the inclination of the rotating 
CO-disk between the nuclei (Tacconi et al. 1999). Under this assumption and with a radial velocity 
difference of 150kms _1 the true distance between the two nuclei can be calculated. With a position 
angle of the CO-disk of 40° and an inclination of 75° the true separation of the two nuclei is 1.4 kpc, 
their orbital velocity is 155 km s -1 with an orbital period of 27 million years. Assuming a simple 
two-body problem with equal mass of 2 x 10 9 M®, this orbit is unstable with the gravitational 
force being A of the centrifugal force. This means that the two nuclei will separate in the future 
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and probably are already past the pericenter. Whether the system is bound can be estimated 
from the escape velocity. We compute an escape velocity of llOkms" 1 , which is of the same order 
as the radial velocity difference of the two nuclei. From the radial velocity difference the system 
seems to be unbound. However, if the nuclei separate again they will see not only the mass of 
the other nucleus but also the halo of the other galaxy. Conversion of orbital angular momentum 
into spin angular momentum in the halo leads to gravitational braking, the galaxies reach their 
apocenter and fall back together again (Barnes & Hernquist 1996; Mihos & Hernquist 1996). It 
seems therefore likely, that NGC 6240 has recently undergone an encounter and is awaiting its next 
encounter. 



7.2. Tidal tails 

The first indication that NGC 6240 might be an interacting galaxy system was its disturbed 
morphology with loops, branches and arms. Toomre & Toomre (1972) showed that such tails can 
be formed by tidal forces in an interaction of two galaxies. Tidal forces during the interaction pull 
stars from the galaxies into regions at a distance of up to several tens of kiloparsec from the center 
of the galaxies. The formation of tidal tails is therefore an efficient way to dissipate orbital angular 
momentum of the interacting galaxies thus allowing the galaxies to come closer to each other and 
finally merge. That the tails in NGC 6240 are dominated by stellar light can be seen from broad 
and narrow band imaging. The tidal tails visible in the R-band image of Armus, Heckman & Miley 
(1990) are not as prominent in their narrow band image of Ha+[N II] which are the most dominant 
emission lines in this wavelength band. The formation of tidal tails is favoured if at least one 
partner in the collision is subject to a prograde encounter. The spin and orbital angular momenta 
of the nuclei suggest a prograde encounter for what is now the northern nucleus. But not only the 
geometry of the collision is important for the formation of tidal tails. Springel & White (1998) 
and Dubinski, Mihos & Hernquist (1999) show that the mass density distribution of the interacting 
galaxies is an important parameter in the formation of tidal tails. Galaxies with more massive 
halos tend to form less sharp and shorter tidal tails. The tidal tails in NGC 6240 are not as "crisp" 
as in NGC 4038/9 ("The Antennae") or in NGC 4676 ("The Mice"), indicating the collision of 
two massive spiral galaxies. Because the parameter space of galaxy interaction models is too large, 
it is impossible to derive a detailed history of NGC 6240 from the morphology of its tidal tails. 
However, simulations of interacting and merging galaxies show, that tidal tails can only be formed 
after the first encounter of the galaxies. Because we see tidal tails in NGC 6240 which are quite 
long, it must have undergone at least one encounter. The tidal tails of NGC 6240 were most likely 
created in the first encounter. How many passages NGC 6240 has undergone we cannot determine 
from the tidal tails. But because the two nuclei are still distinct from each other, it seems that 
only few encounters can have happened and that NGC 6240 is in a rather early merger phase. 

Doyon et al. (1994) showed from K-band photometry that the overall brightness distribution 
of NGC 6240 is very similar to that of an elliptical galaxy. They concluded that NGC 6240 is in 
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an advanced merger state and in the process of forming an elliptical galaxy. We believe that this 
evidence is rather weak and the conclusion should be treated with caution. 

1. The near-infrared light of the nuclei comes from supergiants and not from the old stars 
representing the overall stellar population. Population and L/A4 changes with radius are 
rather likely. 

2. While a r 1 / 4 -law or a King-profile fits the data it is not unique to elliptical galaxies. 

3. The stellar dynamics definitely shows that the two galaxies are still independent entities and 
have not yet merged. 

However, the small separation of the nuclei and the decoupling of the gas from the stars in 
NGC 6240 strongly suggest that the nuclei will merge. But as the formation of tidal tails depends 
on the mass density distribution of the interacting galaxies the merging time also varies with the 
mass density distribution. Galaxies with extended massive halos have longer merging times than 
galaxies with compact low mass halos. Only in the final merging phase violent relaxation leads to 
a strong and rapid dynamical evolution also on global scales (Mihos & Bothun 1998). 

What type of galaxy the merger remnant will resemble we cannot predict from the data. 
However, NGC 6240 might form an elliptical galaxy of which it shows already signatures. The 
stellar velocity dispersion in the nuclei of NGC 6240 is typical for an elliptical galaxy. Outside 
the nuclear region the K-band light can be fit by a r^-law which also approximates the surface 
brightness profile of merger remnants in numerical simulations of galaxy mergers (Barnes 1992). 

7.3. Nuclear Starbursts 

Another hint at how many passages NGC 6240 has already undergone comes from the age and 
duration of the nuclear starbursts. In §4.2 we determined the age of the starburst to be 15 to 25 
million years. Numerical simulations predict that such nuclear starbursts in merging galaxies are 
triggered by close passages of the two galaxies. If this was true, the last encounter happened ~ 20 
million years ago. Because the tidal tails are quite long and pronounced, the first encounter must be 
well past and the last passage was probably at least the second encounter. The dynamical time-scale 
of the orbital motion of the two nuclei is a few ss 10 million years and is of the same order as the age 
of the starburst. The extent of the superwind (« 5kpc) and its velocity 500 km s -1 ) yields an 
age for the superwind of ~ 10 million years, similar to the lifetime of a superwind (Heckman et al. 
1990). If taken into account that the superwind sets in roughly 10 million years after the starburst, 
the superwind and the K-band continuum can be explained by the same starburst. The duration 
of the starburst is much shorter than the age of the starburst and can be explained by a negative 
feedback effects of the starburst through supernovae and stellar winds. Numerical simulation by 
Mihos & Hernquist (1996) predict an episodic star formation during the merging of two galaxies. 
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Gas streaming into the central regions gets compressed and triggers star formation. Depending on 
the existence of massive bulges one major starburst takes place either after the first encounter or 
during the final merging of the two galaxies. Because NGC 6240 is clearly not in the stage of final 
merging, the starburst we see has likely been triggered in an early encounter. 



7.4. Gas Concentration between the Nuclei 

The presence of a self-gravitating cold molecular gas concentration between the nuclei of 
NGC 6240 is not expected from numerical simulations of galaxy mergers (Barnes &; Hernquist 
1996; Mihos & Hernquist 1996). These simulations predict that the highly dissipative gas follows 
the stellar distribution. Gas disks form around the individual galaxy nuclei and get more compact 
as the merger advances. From the point of view of these simulations it would appear as if NGC 6240 
has to be in a rare transient phase. 

However, NGC 6240 is not the only galaxy with a gas concentration between two nuclei. 
VV 114 (Yun, Scoville & Knop 1994), NGC 6090 (Gao et al. 1998; Bryant & Scoville 1999) and 
NGC 4038/9 (Stanford et al. 1990) all have prominent CO-emission peaks between the two nuclei. 
In all cases the projected distance between the nuclei (6kpc for VV 114, 3.5 kpc for NGC 6090 
and 7kpc for NGC 4038/9; Ho = 75kms _1 Mpc -1 ) is much larger than the separation of 750pc 
for NGC 6240. NGC 6240 is also an exception because very little CO-emission is observed at the 
position of the nuclei. In NGC4038/9 three clearly distinct gas concentrations are detected: two at 
the position of the nuclei; and a third between the nuclei, which contains roughly half of the total 
gas mass. NGC 6240, too, has half of the molecular gas in a central concentration between the 
nuclei but the other half is in filaments extending from the central gas concentration. In VV 114 
the molecular gas is in a bar-like concentration with tails extending from it. NGC 6090 shows 
a ridge-like structure of molecular gas which covers also both nuclei. VV 114, NGC 6090 and 
NGC 4038/9 are all, mainly based on the large separations, considered to be early mergers. The 
small separation of the two nuclei in NGC 6240 could be interpreted as if NGC 6240 was in an 
advanced merger state. However, to derive the evolutionary state of merging galaxies only on the 
nuclear separation can be misleading. Due to projection effects and unknown orbits of the galaxies 
the nuclear separation could be much larger. 



8. Conclusions 

From the high resolution near-infrared integral field spectroscopy of NGC 6240 we find: 

• The K-band light of the two nuclei in NGC 6240 is dominated by red supergiants which must 
have been formed in a starburst 15 — 25 million years ago. The duration of the star formation 
is < 5 million years, thus only a small fraction of the starburst age. The total mass of stars 
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formed in the starburst is A4* = 0.4 - 2 x 10 8 M & - 

• The stars in the two K-band peaks exhibit fast rotation. From the stellar velocity field the 
dynamical mass of the nuclei is determined to be M.^ yIi = 2 — 8 x 10 9 M® within the central 
500 pc with similar masses for both nuclei. This exceeds the mass of the most massive star 
forming regions by more than 100, implying that the two infrared emission peaks are the 
massive, rotating bulges of two interacting and merging galaxies. 

• After correction for the rotation component the stellar velocity dispersion peaks between the 
two nuclei at the emission peak of the molecular gas. The inferred mass concentration is a 
self- gravitating concentration of molecular gas. 

• The H2 emission of NGC 6240 peaks between the two nuclei and is thermally excited in a 
slow, continuous shock triggered by the collision of the two galaxies. Molecular gas streams 
to the centers of the galaxies where the turbulent motions also give rise to shock excited H2 
emission. 

• From the orbital and spin angular momenta of the interacting galaxies it seems that one 
galaxy in the system is subject to a prograde encounter. Based on the morphology of the 
tidal tails, the starburst scenario, the kinematics and the excitation of the H2 emission, we 
argue that we observe NGC 6240 shortly after an early encounter which triggered the observed 
starburst. 

The authors are grateful to the staff of the ESO 2.2 m telescope and the Anglo- Australian Tele- 
scope for their support during the observation. We also would like to thank S. Mengel, J. Gallimore, 
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and ROGUE and collecting the data. 
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Fig. 1. — Selected line maps of NGC 6240. The crosses in all four panels denote the position of 
the K-band nuclei. Contours are in 10% steps. Top-left: Integrated continuum emission in the 
wavelength range from 2.2/xm to 2.4/im. Top-right: Map of the CO 2 — absorption bandhead 
depth. Bottom-left: Brackett 7 emission line map. Bottom-right: Integrated H2t> = l — 0S(1) 
emission line map. 
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Fig. 2. — R = 2000 spectra of the H2 v = 1 — S(l) emission peak, the northern nucleus and southern 
nucleus of NGC 6240 (each in an 1" diameter circular aperture). 
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Fig. 3. — Line-of-sight velocity profiles of the H2 v = 1 — S(l) and CO J = 2 — 1 molecular emission 
lines at a velocity resolution of 140 kms -1 . The dashed line represents the H2 v = l — S(l) profile, 
the solid line the CO J = 2 — 1 profile. The aperture size over which the profiles were integrated is 
(y/6 x (y/6. The profiles at the position of the two nuclei are indicated by the labels "North" and 
"South". 
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Fig. 4. — Extinction map of NGC 6240 in contours over a grayscale representation of the K-band 
map. The numbers give Ay for the foreground screen model. 
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Fig. 5. — Spectrum of the southern nucleus of NGC 6240 together with spectra of four giant and 
supergiant stars. For convenient comparison the spectrum of NGC 6240 was blueshifted into the 
restframe of the stellar spectra. The stellar spectra were convolved with a Gaussian profile with a 
FWHM corresponding to 600 km s -1 . 
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Fig. 6. — Results from starburst simulations as a function of starburst age. Shown are four starburst 
scenarios with continuous star formation and starburst durations of 1, 5 and 20 million years. The 
vertical hatched bar denotes the age range from 15 to 20 million years, the age of a red supergiant. 
Top: Brackett 7 equivalent width as a function of starburst age. The measured equivalent widths 
of the southern and northern nucleus are denoted by horizontal hatched bars. Middle: K-band 
luminosity to stellar mass ratio as a function of starburst age. Bottom: Bolometric luminosity to 
K-band luminosity ratio as a function of starburst age. The horizontal hatched bar denotes the 
value for NGC 6240 in a 5" diameter circular aperture. 
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Fig. 7. — Stellar velocity field of NGC 6240 in contours over a gray scale representation of the 
NICMOS K-band map. The contours are spaced 50kms _1 , solid contours represent positive, red- 
shifted velocities, dashed contours represent negative, blueshifted velocities. The systemic velocity 
used is as 7300 km s -1 . 
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Fig. 8. — Measured and model rotation curves for the northern (right panel) and southern (left 
panel) nucleus of NGC 6240. Open diamonds represent the measured rotation velocities, the solid 
line is the best fitting model rotation curve, and the dashed lines denote the la error of the model 
rotation curves. 
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Fig. 9. — Stellar velocity dispersion of NGC 6240 along a north-south direction. The error bars 
along the abscissa represent the width of the aperture, over which the spectra were integrated. 
Velocity error bars are la error bars. 
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Fig. 10. — Channel Maps of the B.2V = 1 — 0S(1) emission line at velocity intervals of 70kms~ 1 . 
Contour levels start at 5 x 10 -6 ergs^ 1 cm -2 st^ 1 in intervals of 5 x ergs -1 cm~ 2 st -1 . In the 
upper left corner of each panel the radial velocity with respect to 7300 km s -1 is plotted. Spatial 
resolution is 0"75 x 0775. 
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Fig. 11. — Schematic view of the assumed collision geometry in NGC6240. For the drawn geometry 
the receding northern nucleus is further away than the southern nucleus. 
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Table 1. Broad and Narrow Band Luminosities of NGC 6240 



Wave Band 


Aperture 


Diameter 


L 






(1) 


(2) 


(3) 


(4) 




(5) 






1 


5.8 x 


10 s 


6.8 x 10 8 




S 


1.5 


1.1 X 


10 9 


1.2 x 10 9 






2 


1.5 x 


10 9 


1.7 x 10 9 


K-Band 
















1 


2.2 x 


10 8 


2.3 x 10 s 




N 


1.5 


4.5 x 


10 8 


4.7 x 10 8 






2 


7.0 x 


10 8 


7.3 x 10 8 


H 2 v = l-0S(1) 


H 2 


4 


2.9 x 


10 7 


3.0 x 10 7 


H 2 f = l-0S(0) 


H 2 


4 


9.3 x 


10 6 


9.8 x 10 7 


H 2 v = 2-lS(l) 


H 2 


4 


4.0 x 


10 6 


4.2 x 10 7 


H 2 f = 2-lS(2) 


H 2 


4 


1.8 x 


10 6 


1.9 x 10 7 




S 


1.5 


9.5 x 


10 5 


1.1 x 10 6 


Brackett 7 


N 


1.5 


8.0 x 


10 5 


8.3 x 10 5 




N3 


1.5 


4.5 x 


10 5 


4.5 x 10 5 



Note. — Columns: (1) Wave band; (2) Position of circular aper- 
ture: southern nucleus (S), northern nucleus (N), H 2 v = l — S(l) 
emission peak (H 2 ) and radio peak N3 (see Colbert et al. 1994); (3) 
Diameter of circular aperture in seconds of arc; (4) Luminosity in 
L© (Ho = 75kms~ 1 Mpc -1 ); (5) Dereddened K-band luminosity 
in Lq. Extinction values are for a foreground screen model from 
53.2. 
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Table 2. Stellar and Dynamical Masses of NGC 6240 



Nucleus 
(1) 


(2) 


i 

(3) 


V ro t 

(4) 


Mdyn 

(5) 




M d y n /M* 

(6) 


North 
South 


0.4 - 1.2 x 10 8 
0.8 - 2.3 x 10 9 


45° 
60° 


360 ± 195 
270 ± 90 


™i£ x 

4.3l|| x 


10 9 
10 9 


> 14 

> 8 



Note. — Columns: (1) Nucleus; (1) Stellar mass in M & within r < 235 pc; 
(2) Stellar mass in A4q within r < 235 pc; (3) Inclination of stellar disks; (4) 
Rotation velocity in kms -1 ; (5) Dynamical mass in Mq within r < 235 pc; 
(6) Ratio of dynamical and stellar mass. 



